Surface modes and critical velocity in trapped Bose-condensates. 
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We propose the novel mechanism leading to superfluidity breakdown in dilute Bose-condensed 
gases. We discuss the properties and highlight the role played by surface excitations in trapped 
condensates. We show that the critical velocity measured in recent experiments is nothing else but 
the critical velocity associated with the spontaneous creation of surface modes. The latter turns 
out to be larger than that found by the analysis of a superflow stability with respect to vortex 
nucleation. 



Superfluidity is one of the most striking manifestations 
of quantum statistics in many-body interacting systems. 
It is usually defined as "the property of a liquid to move 
without friction inside capillaries or around obstacles" . 
Superfluidity is confined to small flow velocities. There 
is a certain critical value of the fluid speed v c beyond 
which the dissipation settles in. The value of the critical 
velocity is directly related to excitations properties . 
The analysis of the energy conservation in a moving liq- 
uid shows that the flow becomes thermodynamically 
unstable when the liquid velocity exceeds 
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Here e(p) is the energy of an excitation with the momen- 
tum p in the laboratory frame. 

When applied to liquid helium, the criterion (jlj) gives 
the critical velocity coinciding with the sound velocity, 
which is too high to explain the experimentally measured 
value. Moreover, the experimental results appear to be 
non-universal and depend on the geometry of the liquid 
sample J| . This contradiction was resolved by R. Feyn- 
man who suggested, that the dissipation is related 
to creation of quantized vortices. If this is taken into 
account, the critical velocity turns out to be both suffi- 
ciently small and depends on the external factors. The 
combination of the vortex hypothesis and the Landau 
criterion ([!]) very well explains the experiments in liquid 
helium j|. 

These fundamental theories are now being tested in the 
new experiments with ultra-cold Bose-condensed gases 
(see 0|| for a review of experimental and theoretical ac- 
tivity in the field of Bose-Einstein condensation (BEC) 
in weakly interacting gases). It is now possible to study 
dynamical properties of superfluidity in systems, which 
are both exceptionally well controlled experimentally and 
understood theoretically. The closest analogy for the 
critical velocity measurement in long capillaries is the 
experiment done at MIT Q . There the onset of the dis- 
sipation caused by a moving laser beam was studied as a 
function of the beam velocity. The reported value of the 
critical velocity is well below the speed of sound. This 
can be explained either by the vortex nucleation |||7| , or 
by the modification of the excitation spectrum in inho- 
mogeneous condensates H. 



Another way to measure critical velocity is to rotate 
a superfluid and find the minimum (critical) angular ve- 
locity f2 c at which an excitation shows up in the system. 
The straightforward generalization of Eq.(^) for the case 
of a rotating superfluid gives : 
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where e(M) is the energy of an excitation as a function of 
its angular momentum M. In a Bose-condensed gas, the 
lowest value of f2 c is obtained when e(M) is the energy 
of a single vortex state. Remarkably, the critical velocity 
(||) is smaller than that observed at ENS ||. 

In this Letter we perform a simple analysis of criti- 
cal velocity issues in a trapped Bose-condensed gas. We 
point out that the universal criterion (Q) gives no infor- 
mation about the dissipation time scales. Actually, Eqs. 
(|]) and (||) merely reflect the fact that dissipation is not 
forbidden by the energy conservation law as soon as the 
(angular) velocity of the liquid reaches a certain value. In 
spite of the fact that vortex excitations may have the low- 
est values of the critical velocity, in a dilute gas they can 
hardly be formed suddenly out of nothing. Instead, as it 
is prompted both by the experiment [ fLCfl and numerical 
simulations | pT| , vortices are nucleated at the conden- 
sate border from the condensate shape oscillations. This 
can only happen if the velocity of the superflow exceeds 
the critical velocity corresponding to the formation of 
the necessary surface modes, which later evolve into the 
vortices in the course of non- linear condensate dynamics. 
To support this view we analyze the elementary excita- 
tions of a confined Bose-condensed gas. We show that, in 
a sufficiently dense condensate, the lowest possible val- 
ues of the critical velocity correspond to the excitations 
localized at the low density region close to the edge of 
the condensate. Remarkably, the surface modes turn out 
to be responsible for the previously reported decrease of 
v c compared to the velocity of sound in inhomogeneous 
condensates We note that this effect is a novel fea- 
ture characteristic of dilute Bose-condensed gases and is 
absent in "conventional" superfluids, such as liquid he- 
lium. In the latter case the superfluid is incompressible 
and the surface excitations of the mentioned type cannot 
exist. Using a simple model of the surface modes, we cal- 
culate both v c (for an MIT-type experiment) and Q c (for 
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an ENS-type experiment), assuming that the generation 
of the surface modes is the necessary precondition for the 
vortex nucleation. 

Consider a Bose-condensate confined in an infinitely 
long cylindrical harmonic trap characterized by the trans- 
verse frequency to (hereafter we assume Ti — m = 1, where 
in is the mass of the condensate atom) . The condensate 
is characterized by its chemical potential u and is as- 
sumed to be in the Thomas-Fermi(TF) regime: u ^> to. 
In this case, the condensate density has a simple form: 
n o(p) = "o(l — P 2 /Rc)i where p is the radial coordinate 
and R c = (2^/w 2 ) 1 / 2 is the TF size of the condensate. 
The condensate density on the axis of the trap no is re- 
lated to the chemical potential: p = 4-7moa, where a is 
the scattering length || . 

The character of the excitations in this geometry is 
pretty well understood by now. The lowest elementary 
excitations are phonons. The velocity of sound is found to 
be \[2 times smaller than that in a homogeneous conden- 
sate of the same density n Q : cz = y/gn a /2 pl|. The ra- 
dial profile of the phonon modes turns out to be very sim- 
ilar to the profile of the condensate wavefunction. Hence, 
the phonons are delocalized inside the condensate and are 
analogous to the "bulk" phonons in homogeneous Bose- 
condensates or in liquid helium. As the axial wave- vector 
of the excitations p increases (i.e. when pR c ~ 1), the 
radial dependence of the phonon modes becomes more 
and more involved, and the spectrum ceases to be linear 
|l2"| . The behavior of the excitations for arbitrary values 
of p was studied numerically in . The critical velocity 
obtained with the help of Eq. (jl|) was shown to decrease 
together with the TF parameter rj = ui/fj, and can be 
much smaller than the velocity of sound cz- 

As we will see below, the excitations with the lowest 
critical velocity are in fact the surface modes, i.e. the 
excitations localized at the TF border of the condensate. 
In this region the condensate wavefunction ipo has a uni- 
versal form and can be analyzed with the help of the 
following universal equation ||l3| , [L4[ : 
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Here x = (p - R c )/d, d = i? c (w 2 /2u 2 ) 1 /3 ^ r c j s 
the "condensate border thickness", and xo — {d 2 g) 1 ^ 2 ipo 
is the dimensionless condensate wavefunction. Deeply 
inside the condensate spatial region ( x — > —00), the 
solution xo matches the TF density profile: xo 
yf—x. Far away from the condensate border [x —> 
00 ) the condensate density quickly decreases: Xo ~ * 
0.19exp(-2a; 3 / 2 /3)/x 1 / 2 f|. 

The spectrum of elementary excitations can be found 
by studying the behavior of quantized linear fluctuations 
of the order parameter Sip — if} — tpQ, where ip is the 
annihilation operator of the atomic field Q , and 
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Here d pm ,{d) prn ) are the annihilation(creation) operators 
of the Bogolyubov excitations characterized by the lon- 
gitudinal and angular momentums p and m, respectively. 
One can show that the functions f^ m = u pm =F v pm satisfy 
the following Bogolyubov-de Gennes equations [^5[: 
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The whole spectrum of the elementary excitations can 
be recovered by solving Eqs.([I]),(|^) numerically together 
with the exact Gross-Pitaevskii equation for the conden- 
sate wave function ipQ. Below we will take advantage of 
the representation (0) and analyze only the modes con- 
fined to the vicinity of the condensate border \p—R c \ ~ d. 
This can be done by rewriting Eqs.(^),(|^) in the same 
units as Eq.(^). Since the excitations are radially lo- 
calized inside a shell with the thickness d <C R c , the 
Laplacian A can be substituted by the second derivative 
in x, and the centrifugal term can be approximated by 
its value at p = R c . Thus, Eqs.@,@ take the simple 
universal form 
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where e = ed 2 , and 
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Since Eqs.(||),(0) contain only the single parameter fc 2 , 
the spectrum of the excitations is given by the universal 
relation: e pm = d~ 2 F(k 2 ). Here F is a dimensionless 
spectral function found numerically (see Fig. nl). 
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FIG. 1. Dimensionless spectral function F for the lowest 
surface excitations obtained by numerical solution of Eqs.(^j) 
and (|). 
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To make sense of the numerical results we resort to a 
simple model, based on WKB approximation to Eqs.(||) 
and (Q). Quasiclassical radial momentum p p has the fol- 
lowing Bogolyubov-type form 
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To make the things even more simple, we ignore the ex- 
act behavior of the condensate wave function close to 
the condensate border |x| % 1 and substitute xo with 
its asymptotic values at \x\ ^> 1. For sufficiently high 
momentums (k 2 ^> 1), we find, that the classical turning 
points are close to the condensate border x = 0. The di- 
mensionless energy of the excitations is much larger than 
the mean field interaction Xo- Hence, the dispersion re- 
lation rti) can be written as 
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Born-Sommerfeld quantization rule gives essentially the 
single-particle spectrum 
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where n = 0,1,... is the radial quantum number, and 
C n = (37r/8\/2) 2 / 3 (4n + l) 2 / 3 is the dimensionless con- 
stant found from the boundary condition at the clas- 
sical turning points |l(||l7j]. The quasiclassical value 
Co = 0.885 found in our model calculation is remark- 
ably close to that obtained by exact numerical solution of 
Eqs.(||) and (^): Co ~ 0.9. This coincidence is somewhat 
surprising, since our analysis fails to recognize the com- 
plicated character of the condensate wavefunction close 
to the condensate border. At the same time the numeri- 
cal solution shows that the mode functions are localized 
at |x| ^ 1, where our approximations are better justified. 
Moreover, WKB approximation is known to provide very 
reasonable estimations for excitation energies even if the 
quantum numbers are small [[l7f . The exact solution of 
Eqs. (Q) and (0) also shows that the approximation (||) 
holds fairly well even for pd,md/R c ~ 1, in spite of the 
fact that the excitations characterized by small momen- 
tums k < 1 spread more inside the condensate spatial 
region and, hence, are not single particles (see Fig. |l|). 

We note that, although the surface modes (||) arise 
from the analysis of the Bogolyubov dispersion relation, 
they are quite different from the Bogolyubov excitations 
in a homogeneous condensate. At smaller momentums 
(pd < 1) the localization size of the solutions to Eqs. 
@ > (0) increases and gradually approaches the size of the 
condensate. In this limit, the dispersion relation recovers 
its phonon limit. The surface modes (^|) exist in a thin 
layer around p = R c and which is still thicker than the 
healing length Iq = l/^/8nan~o~ in the condensate: rf> Iq 
in the TF limit. This may happen only in a compress- 
ible liquid and thus the dispersion relation similar to (||) 



can hardly be found in "traditional" superfluids, such as 
liquid helium. In the latter case the condensate density 
drops on a distance scale of order Iq at the wall of the 
confining vessel. 

The validity of Eq. (|]) is restricted by a number of as- 
sumptions. First, we require that the radial localization 
length of the surface excitations be much smaller than 
the size of the condensate. Since we are mostly inter- 
ested in e pm ~ 1/rf 2 , the radial size of the modes is 
~ d <C R c in the TF regime. In the same way we can 
neglect the spatial variations of the centrifugal potential 
provided that the momentums of the involved excitations 
are sufficiently small. Indeed, a simple WKB calculation 
shows that the first correction to Eq. (JsJ) arising from 
the variations of the centrifugal potential close to p w R c 
results in a slight change of the ground state energy: 
ACo/Co ~ m 2 p/R 2 . The solution of Eq.(||) corresponds 
to m ~ R c /d <C Rc/p 1 / 2 , so that for our purposes, the 
finite-m corrections to C n can be safely neglected. 

The dispersion relation (||) is the central result of our 
work. It allows us to study the critical velocities (|l|) and 
(0) in the deep TF regime (rj <C 1). First we turn to 
an MIT- type experiment. The application of criterion 
([!]) shows that the lowest value of critical velocity corre- 
sponds to the instability of the excitations with m = 0: 
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This result has to be compared with the stability analy- 
sis of the vortex excitations. Consider a superfluid flow 
of velocity v s without vortices. The energy of the liq- 
uid scales as ~ p s v 2 R 2 L, where p s no is the super- 
fluid density and L is the length of the trap. Let us 
consider instead the vortex configuration similar to that 
discussed in relation to the experiment Jl8[ , i.e. a se- 
quence of vortices at the distance tt/v s from each other 
along the trap axis. The average velocity of the super- 
fluid would be the same v s . The energy of this vortex 
chain is roughly the sum of the vortex energies and thus 
can be estimated as p s R c \og(\/ri)v s L. Comparing the 
energy of the superflow with and without the vortices, 
we find that the vortex configuration has a lower energy 

as soon as v s > ~ r/log(ri~ 1 )cz ■ This means that 
in a TF condensate, vortices form the true ground state 
already at v > Vc V \ At the same time, since the vor- 
tices are nucleated only at the edge of the condensate, 
they originate from the surface excitations, which only 
become unstable if v > v c > ■ Therefore, in the TF 
regime no vortices can be nucleated unless the superfluid 
velocity reaches v c . This statement does not contradict 
the conclusions derived in |Q, since the vortex states 
posses very small energies compared to the irrotational 
superflow of velocity v c , and therefore further relaxation 
leads to the creation of vortices. 

The dispersion relation for the surface excitations also 
gives us a way to discuss the critical angular velocity 
in an ENS-type rotating superfluid experiment. Again, 
applying the criterion (0) to the dispersion relation (m we 
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immediately find that the surface modes in a cylindrical 
trap first become unstable if p = and 
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This value has to be compared with the lowest angular 
velocity at which a vortex can be formed in the conden- 
sate. It can be derived using the same kind of simple 
arguments as above: fic ~ ujrjloglr/^ 1 ) ffl. As in the 
discussion above, in the deep TF limit Qc < £l c and, 
we conclude that, in order to make up a vortex on the 
surface of a rotating superfluid, one has to first reach the 
angular velocity ([n]). Only then the surface modes are 
nucleated and, finally, a vortex appears in the superfluid 
in the course of the turbulent liquid dynamics. This point 
of view can be supported by earlier work JlO| , where the 
authors conclude that f2 c coincides with the frequencies 
of the phonon-like excitations with m = 2, 3 (depend- 
ing on the symmetry of the rotating perturbation). We 
argue that the TF parameter in the ENS experiment is 
not sufficiently small. Since the critical velocities jltj| ) 
and ([ll]) change very slowly as rj decreases, the distinc- 
tion between the bulk and the surface modes is not very 
sharp at moderate values of the TF parameter. 

Finally, we would like to comment on the physical na- 
ture of criteria ^ and (||). Consider an obstacle (probe) 
of a size I moving inside a superfluid with the velocity v. 
In a weakly interacting Bosc-condensed gas the dissipa- 
tion arises pair creation of excitations with radial momen- 
tums p and p' . The energy conservation law |Q requires 
that 
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If the fluid velocity slightly exceeds v c (Q), both solutions 



of Eq. (|12j) p, p' where p* is found from the equation 

e(p*) = Vcp*. At the same time, the uncertainty principle 
requires (p +p') 1/2- Therefore, pj, < 1 and thus the 
dissipation at v « v c is only possible if the probe is suffi- 
ciently small: / < P* 1 ■ In our case we find that p* ~ d^ 1 
and thus the probe dissipates only if I < d. In this sense 
the criteria ([!]) and @ give the values of the critical ve- 
locities measured by an ideal, point size probe. We note 
that this requirement is not at all automatically satisfied 
in the experiments. This, in turn, may lead to addi- 
tional complications in ascribing a particular mechanism 
to the reported critical velocity values. For example, a 
finite size probe may nucleate a vortex ring on its sur- 
face. This instability leads to a sufficiently small value of 
the critical velocity: v c ~ l/l (see fjj and refs. therein). 
In the latter case, the mechanism of the critical veloc- 
ity decrease is somewhat similar to that reported here: 
the increase of the superfluid velocity around a moving 
body leads to the decrease of the condensate density, so 
that the moving body-superfluid interface can host an 
unstable "surface" excitation. 

In conclusion, we develop a simple theory of critical ve- 
locity in trapped BEC. We highlight the role played by 



the surface excitations in large TF condensates. We sug- 
gest that, since the vortex excitations are formed at the 
condensate border and thus are initially surface excita- 
tions, the critical velocity measured in the recent experi- 
ments cannot be smaller than the critical velocity associ- 
ated with the generation of the surface excitations. The 
latter turns out to be larger than that found by the anal- 
ysis of thermodynamic stability of superflow with respect 
to vortex nucleation both in a long trap with axially mov- 
ing superfluid and in a rotating trapped Bose-condensate. 
We note that this novel feature is a direct consequence 
of compressibility of trapped Bose-condensed gases and 
hence cannot be related to traditional superfluids, such 
as liquid helium. 
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